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I In the background of the strong magnetic field the vacuum is suggested to possess an electro- 

magnetically superconducting phase characterised by the emergence of inhomogeneous quark- 
antiquark vector condensates which carry quantum numbers of the charged p mesons. The p- 
meson condensates are inhomogeneous due to the presence of the stringlike defects (the p vor- 
tices) which are parallel to the magnetic field (the superconducting vacuum phase is similar to 
the mixed Abrikosov phase of a type-II superconductor). In agreement with these expectations, 
we have observed the presence of the p vortices in numerical simulations of the vacuum of the 
s O quenched two-color lattice QCD in strong magnetic field background. We have found that in the 

quenched QCD the p vortices form a liquid. The transition between the usual (insulator) phase at 
low B and the superconducting vortex liquid phase at high B turns out to be very smooth, at least 
in the quenched QCD. 
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1. Introduction 

Recently it was suggested that in the magnetic field background the vacuum becomes electro- 
magnetically superconducting if the strength of the magnetic field exceeds the critical value [1,2]: 

B c = m 2 p /e^ 10 16 Tesla, (1.1) 

where m p = 775.5 MeV is the mass of the p meson. The magnetic fields which are two-three times 
stronger than the critical value (1.1) are expected to be experimentally reachable at the LHC [3]. 

If the strength B = \B\ of the uniform static magnetic field B = (0,0,5) is higher then the 
critical value (1.1), then the usual vacuum ground state should experience a tachyonic instability 
towards the emergence of a new ground state characterised by the presence of the following quark- 
antiquark condensates: 

(und) = p(x ± ) , (uy 2 d) = -ip(x ± ) . (1.2) 

The complex scalar field p is a function of the (transverse) spatial coordinates x± = (x\,X2). The 
quark- antiquark condensates (1.2) carry the quantum numbers of the electrically charged p mesons 
so that this phenomenon may also be interpreted as the p -meson condensation. 

Since the condensed pairs (1.2) are electrically charged states, their condensation implies, 
almost automatically, the emergence of the electromagnetic superconductivity in the new vacuum 
state at B > B c . The emerging superconducting vacuum state has quite unusual properties (for a 
detailed review, see Ref. [5]): 

(i) Anisotropy: the ground state is a perfect conductor for the electric currents directed strictly 
along the magnetic field axis. In the transverse directions the superconductivity is absent. 

(ii) Inhomogeneity: the transport coefficients depend on the transverse coordinates x±. 

(iii) Absence of the Meissner effect : the induced superconductivity cannot screen the background 
magnetic field due to the mentioned spatial anisotropy. 

(iv) In-tandem superfluidity: one can argue that the superconducting state (1.2) should always be 
accompanied by a superfluidity of the neutral p(°) mesons [1,4]. 

(iv) Optically, the vacuum superconductor is metamaterial with a perfect lens properties [6]. 

In the magnetic-field-induced vacuum superconductivity, the quark-antiquark composites (1.2) 
play the same role as the Cooper pairs play in a conventional superconductor. This analogy may go 
even deeper: in a very strong magnetic field a conventional type-II superconductor was suggested to 
enter a quantum limit of "reentrant superconductivity" characterised by a p-wave spin-triplet pair- 
ing, absence of the Meissner effect, and a superconducting flow along the magnetic field axis [7]. 
It is encouraging that these are exactly the features which we expect to be realised in the supercon- 
ducting vacuum state at B > B c . 

There are indications from holographic [8, 9] and numerical [10] approaches that the p-mesons 
should be condensed in the strong magnetic field; see also the ongoing discussion in Refs. [11, 12]. 
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2. Vortex lattice ground state in the mean-field approximation 

The p -meson condensation can be characterised by a single scalar function, 

p(x) = ^u(x)y + d(x) , 7+ = 71 + iy 2 , (2.1) 

where the combination y + of the Dirac matrices corresponds to the s z = +1 projection of the p- 
meson spin onto the magnetic field axis. In the small-condensate limit, |p| <C m p , the ground state 
of the p -meson condensates can be described by the following general form [1, 2, 4, 9]: 

{jr 7 ^ 'ZTC 

- W72 ( kl 2 + z 2 ) - ^ v v + 2nvn -r r ' Lb = v -s ' (2 - 2) 

where L# is the magnetic length and V is a real parameter. The solution (2.2) is similar to - and 
inspired by - the Abrikosov vortex lattice configuration which appears in a mixed state of a type-II 
superconductor in the magnetic field background [18]. 

In Eq. (2.2) the complex coefficients C n are fixed by the energy minimisation condition. It is 
usually assumed that the coefficients C n obey the N-fo\A symmetry, where N is an integer [18]: 

C n+N = C n , N =1,2,... (2.3) 

so that the condensate (2.2) has a periodic structure in the transverse (jq ^2) plane. A condensate of 
the form (2.2) possesses an infinite set of zeros, which mark the centres of the topological string de- 
fects called the p vortices [1]. Similarly to the Abrikosov vortices in conventional superconductors, 
the condensate p acquires the phase shift 2n as we circumvent these zeros. 




ordered (lattice) vortex state perturbed (liquid/gas) state 



Figure 1: (left) The equilateral triangular vortex lattice in the mean-field approximation to the ground state 
of QCD and (right) its suggested melting due to the presence of quantum and/or thermal fluctuations. 

The simplest configuration with N = 1 (all C w 's are equal) corresponds to the square vortex 
lattice. However, the true mean-field ground state is given by an equilateral triangular vortex lattice 
(which is sometimes also called "hexagonal lattice") with the following set of parameters [4, 9, 18]: 

N = 2, Ci=±iCb, v = ^3/^2^0.9306. (2.4) 

Thus, in the mean field approximation, the p -meson vortex state is similar to the true ground state 
of the Abrikosov lattice in a type-II superconductor [18]; see Fig. 1 (left) for an illustration. 
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The analytic results of Ref. [4] suggest that the p vortices are weakly interacting with each 
other in the vortex lattice state. This fact means that the presence of thermal and/or quantum 
fluctuations may either melt the vortex lattice state to a liquid vortex state or even evaporate it by 
forming a vortex gas (such phenomena are known to happen with the Abrikosov vortex lattices in 
conventional superconductors [14]). The vortex liquid state is a superconducting state characterised 
by the broken t/(l) e .m. electromagnetic gauge symmetry while the vortex gas state is, generally, 
a normal (non-superconducting) state with the unbroken t/(l) e .m. group. In the context of the 
solid state physics, the vortex lattice-liquid-gas phase diagram for magnetic-field-induced reentrant 
superconductivity was discussed in Ref. [13]. 

3. Numerical simulations 

In our numerical setup we basically follow Ref. [10]. We use lattice Monte-Carlo simulations 
of SU (2) Yang-Mills lattice gauge theory. The quark fields are introduced by the overlap lattice 
Dirac operator <2) with exact chiral symmetry [16], and - due to the presence of the magnetic 
field - with twisted spatial boundary conditions [17]. The quarks are treated in the quenched 
approximation so that the vacuum quark loops are absent in our approach. We have studied 18 4 
and 19 4 lattices in a wide range of the magnetic field strengths, eB = (0 ... 2. 14) GeV 2 with lattice 
spacings a « 0. 1 1 fm. We have used 20 gauge configurations per each value of the magnetic field. 

An explicit manifestation of the superconducting phase would be the presence of the conden- 
sate (1.2). Unfortunately, the observable (1.2) cannot be computed directly in our approach, while 
it can be accessed via the following simplest p -meson correlator: 

♦ W ^ - (P'(O)PM )„„ ^ Tr ( %{A l B) +m y, %(A l B) — ») , (3.1) 

where the p meson field is defined by Eq. (2.1). The subscripts in Eq. (3.1) indicate that the 
correlation function (jc) is computed in the fixed background of both the non- Abelian gauge field 
A and the Abelian magnetic field B. Equation (3.1) represents this vector correlator in terms of the 
(overlap) Dirac propagators in the background of both Abelian and non- Abelian gauge fields. 

We notice that under the Abelian transformation from the electromagnetic gauge group £/(l) e . m 
the field (3.1) transforms as a charged scalar field 1 : <j>(x) — > e ieco ^(j)(x). However, the effective field 
(x) is still a two-point correlation function which falls off exponentially as the distance x increases. 
This property is not a desired behaviour for a genuine local scalar field so that the quantity (3.1) 
cannot, strictly speaking, be associated with the p -meson field itself. 

Fortunately, we may get an insight from Ref. [15] where a qualitatively similar issue was 
encountered. In that work the chromoelectric flux tube was studied using a rectangular Wilson loop 
W as a source and the local energy density operator & as a probe. Although the expectation value 
of the Wilson loop falls off exponentially as the area of the Wilson loop grows, the energy density 
in the presence of the Wilson loop, given by the normalised energy ratio {^)^ = ( 0^ } / (W), 
has, generally, a non- vanishing profile as the area of the Wilson loop grows. 

The gauge transformation at the origin, (p(x)^e ie0) W(p(x), acts as a global phase which is not essential for our 
interpretation of the effective field <f>(x). 
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By analogy with Ref. [15], we consider the normalised scalar energy of the p -meson field 
E(x), the normalised electric (super)current j^(x) generated by the p -meson field, and the local 
vortex density v(x) in transversal (x,y) plane, respectively (we use the continuum notations in 
order to simplify the expressions): 

|<H*) 

2i\(j>(x) 



£ W = uL\|2 > Dn = d»- ieA^ , (3.2) 



, N e ab d d 

v(x) smgarg0(x) = — — — arg0(x), cub = 1,2. (3.4) 

l7t OX a OX}y 




r><\- 



(x,y) plane (jc,z) plane 

Figure 2: Two types of the "probe" cross-sections of the (expected) melted p-vortex lattice, Fig. 1 (right). 

In search of signatures of the (perhaps, melted) p-vortex lattice we have studied (configuration- 
by-configuration) the behaviour of the normalised energy density (3.2) in the (x,y) and (jc,z) planes 
(we remind that the magnetic field is directed along the z axis); see Fig. 2. In the center of a 
physical p vortex the energy density is higher than the energy density outside the vortex. Thus, 
if the physical p vortices are formed in the (sufficiently strong) magnetic field background, than 
we may expect the formation of the pointlike lumps of the energy density in the (x,y) plane [see 
Fig. 2 (left)] and the formation of the linelike structures in the (jc,z) plane [see Fig. 2 (right)]. 

Typical examples of the behaviour of the energy density in the (x,y) and (jc,z) planes are shown 
for weak (eB = 0.356 GeV 2 ), moderate (eB = 1.07GeV 2 ) and high (eB = 2.14GeV 2 ) magnetic 
fields in Fig. 3. In accordance with our qualitative expectations, at low magnetic field the vortex 
lattice is not formed. At the moderate magnetic field the formation of a coherent vortex structure 
is seen while the vortices are not strictly ordered in the transversal plane and they are not quite 
parallel to the magnetic field. At higher magnetic field the physical picture is visually consistent 
with the presence of a melted lattice (liquid) of the p vortices; see Fig. 1 (right). 

The peaks in the energy density (3.2) are correlated with the p vortex positions (3.4) and that 
the p vortices are encircled by the supercurrents (3.3). The latter feature is shown in Fig. 4. Thus, 
the numerically observed vortices do indeed bear the essential features of the physical vortices. 

The nature of the p-vortex state may be characterised by the normalised vortex- vortex cor- 
relation function (u(0)u(/?))/(u(0)) , where the p vortex density is given in Eq. (3.4). At 
low magnetic fields this function is a monotonically rising function of the inter- vortex distance R, 
Fig. 5 (left), implying that the p vortices constitute a (nonsuperconducting) gas. 



5 



Vortex liquid in superconducting vacuum of quenched lattice QCD 



M. N. Chernodub 



eB = 0.356 GeV 2 Fluxes = 6 eB = 0.356 GeV 2 Fluxes = 6 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

x, fm x, fm 



eB = 1 .07 GeV 2 Fluxes = 1 8 eB = 1 .07 GeV 2 Fluxes = 1 8 




0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 

x, fm x, fm 



eB = 2.14 GeV 2 Fluxes = 36 eB = 2.14 GeV 2 Fluxes = 36 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

x, fm x, fm 



Figure 3: Typical behaviour of the energy density in the (x,y) planes (the left column) and in the (x,z) 
planes (the right column) for eB = 0.356 GeV 2 (the upper panel), eB = 1.07 GeV 2 (the middle panel) and 
eB = 2. 14 GeV 2 (the lower panel). The number of elementary fluxes is n = eBL x L y / {lit) . In the left column 
the magnetic field is perpendicular to the page and in the right column the magnetic field is directed vertically. 
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0.0 0.5 1.0 1.5 2.0 0.0 0,5 1.0 1.5 

x, fin x, fm 



Figure 4: Examples of the superconducting currents (the blue lines), Eq. (1.6), around the vortices (the red 
squares), Eq. (2.7), in the (x,y) planes at eB = 1.07 GeV 2 (left) and eB = 2.14GeV 2 (right). According to 
the analytical expectations [1,4] the currents should encircle the p vortices in the clockwise direction (an 
example of the mean-field solution, Eqs. (1.4), (1.6), is shown in the inset of the right figure, from Ref. [4]). 



At higher magnetic fields a non-monotonic behaviour manifests itself via the formation of a 
wide maximum at intermediate distances; see Fig. 5 (right). The appearance of the peak indicates 
the presence of the anticipated superconducting p -vortex liquid. 




Figure 5: The normalised vortex-vortex correlation function ( v(0)v(R) ) / ( v(0) } 2 in the (x, y) plane. The 
monotonic (non monotonic) behaviour of the correlator signals the presence of the gas (liquid) vortex state 
at low (high) values of the magnetic field B as shown in the left (right) plot. 

The melting of the vortex lattice in quenched QCD may make it difficult to observe the sug- 
gested p vortex condensation using the standard numerical tools. Indeed, in the vortex lattice 
state [described by Eqs. (2.2) and Eq. (2.4)] the phase of the p-meson field changes by 2k around 
each vortex so that the p -field is an oscillatory function of the transverse coordinates x and y. 
Thus, in the superconducting state at B > B c , the space-averaged (bulk) condensate is always zero, 
(p (*)) bulk = 0' despite the fact that the local p-meson condensate is large and the ground state is a 
superconductor (the same is true for the Abrikosov mixed state in a type-II superconductor [18]). 

However, if the vortices were strictly straight, then the p-meson condensation could in prin- 
ciple still be determined by studying a long-distance limit (taken along the straight vortex world- 
sheets) of the correlation function (3.1) averaged over gluon fields. In this straight- vortex case 
the oscillating phase would not contribute to the correlation function so that long-distance correla- 
tor should generally be nonzero, lim (p(O)p(x")) ~ | (p(0)) | 2 with = (0,0, x,t). However, 
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in the liquid vortex phase the vortex worldsheets are not flat surfaces so that the vortex wob- 
bling may, in general, add large phase fluctuations to the correlator of the p -meson field, hence 
lim (p(O)p(x")) = in this superconducting state. Thus, one may encounter a technical difficulty 

in determination of the exact p -meson condensate in the liquid state by using the p -field correlators. 
4. Conclusions 

We have numerically observed the formation of the p -vortex liquid in the vacuum of quenched 
two-color QCD in strong magnetic field background. The vortex liquid phase is an electromagnet- 
ically superconducting phase characterised by the inhomogeneous order parameter (the p meson 
condensate), similarly to the mixed (Abrikosov) phase of an ordinary type-II superconductor. We 
argue that in this phase the calculation of the (highly inhomogeneous) p -meson condensate by us- 
ing the standard methods should be taken with care. The transition between the usual (insulator) 
phase at low B and the superconducting vortex liquid phase at high B turns out to be very smooth. 
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